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The Hall response provides an important characterization of strongly correlated phases of matter. 
We study the Hall conductivity of interacting bosons on a lattice subjected to a magnetic field. 
We show that for any density or interaction strength, the Hall conductivity is characterized by a 
single integer. We find that the phase diagram is intersected by topological transitions between 
different integer values. These transitions lead to surprising effects, including sign reversal of the 
Hall conductivity and extensive regions in the phase diagram where it acquires a negative sign. 
This implies that flux flow is reversed in these regions - vortices there flow upstream. Our finding 
have immediate applications to a wide range of phenomena in condensed matter physics, which are 
effectively described in terms of lattice bosons. 



The Hall response is a key theoretical and experimen- 
tal tool for characterizing emergent charge carriers [1] in 
strongly correlated systems, ranging from high tempera- 
ture superconductors [3-5] to the quantum Hall effect [2] . 
In this paper, we study the Hall conductivity of strongly 
correlated bosons on a lattice. We find that the entire 
phase diagram of such systems can be characterized us- 
ing a single integer p, and inevitably contains topological 
transitions between different p - values. These observa- 
tions allow us to calculate the Hall conductivity through- 
out the whole phase diagram, and we show that they lead 
to surprising consequences, such as sign reversals of the 
Hall conductivity. The model we study describes a wide 
range of systems in condensed matter physics, to which 
our results have immediate implications. Examples are 
cold atoms on optical lattices [6, 7], Josephson junction 
arrays [8], granular superconductors [9, 10], and perhaps 
even high temperature superconductors such as the un- 
derdoped cuprates [11-13]. 

In the absence of disorder and at weak magnetic fields 
the Hall conductivity of bosonic systems is dominated by 
the flow of superfluid vortices. For a continuum (Galilean 
invariant) superfluid, vortex flow gives a Hall conductiv- 
ity which is proportional to the ratio of the particle den- 
sity and the applied magnetic field. We find that on the 
lattice, vortex dynamics is strongly modified. As a result 
the Hall conductivity is characterized, in addition to the 
particle density, by the integer p. We show how emer- 
gent particle-hole symmetry points in the ground-state 
phase diagram necessarily lead to a non-trivial behavior 
of this integer and we discuss the topological transitions 
between different p-scctors. As we shall show, these tran- 
sitions are attributed to degeneracies in the many body 
spectrum, which serve as sources for the Berry curvature. 

Specifically, we focus on the conventional Bose- 
Hubbard model [14] in two dimensions. We restrict our 
study to a dissipation-less system, at zero temperature 
and without disorder. Within the phase diagram of this 
model we find large parameter regions corresponding to 
a negative Hall conductivity, a xy < 0, and reversed vor- 
tex motion where vortices flow upstream, cf. Fig 1. We 
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FIG. 1. Topological transitions in the Bose-Hubbard 
phase diagram. The Galilean invariant regime denotes the 
region where o xy is proportional to the particle density nb 
divided by the magnetic field strength B. Mott insulator lobes 
are indicated in gray. The yellow and green lines exhibit an 
emergent particle hole symmetry, where a xy — 0. They are 
divided into two types: (i) Lines emanating from the tip of 
the Mott lobe at integer boson filling where cr xy has a smooth 
zero crossing (green), (ii) Transition lines (yellow), through 
which Bo xy exhibits integer jumps. The latter continue into 
the phase diagram, with o xy > 0, as indicated by the dashed 
lines. The blue region corresponds to regions where the Hall 
conductivity is negative. 



discuss methods to directly test these results in cold 
atom systems where the neutral atoms are subjected to 
synthetic magnetic fields introduced through rotation or 
phase imprinting [15, 16]. 
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FIG. 2. Forces acting on a vortex, (a) The classical 
Magnus force due to the interaction of the velocity field of 
the vortex and the external flow v s acts perpendicular to v s . 
(b) Vortex motion leads to a change in the momentum of the 
system due to its phase singularity, which is perpendicular 
to its velocity v„. (c) Moving a vortex around a lattice site 
yields a Berry phase of 2-ira — 2n(nb + p)- 



I. HALL CONDUCTIVITY AND VORTEX 
MOTION 

We begin by giving a description of vortex dynamics 
in bosonic systems. A vortex moving with respect to a 
current experiences a force arising from the interaction 
of the velocity field of the vortex with the one of external 
current. This hydrodynamical force is called the Magnus 
force, and it acts perpendicularly to the current, as de- 
picted in Fig. 2(a). Similarly, a superfluid vortex (of unit 
vorticity) in two dimensions experiences a force 



F M = — 27rfi,n s v s x e 2 . 



(1) 



where n s is the number density of superfluid bosons, and 
v s is their velocity. The unit vector e z is a normal to the 
plane. 

The force F M in Eq. (f) arises from the dynamical 
phase (time integral of the energy) in a Lagrangian de- 
scribing the superfluid. Such a Lagrangian necessarily 
contains a term corresponding to the Berry phase picked 
up by the vortex motion. The Berry phase acquired by a 
vortex moving around a loop of area S is given by 2 / naS 1 
where a is a proportionality factor which depends on the 
microscopic details of the Hamiltonian. Therefore, an 
equation of motion for the vortex leading to dissipation- 
less flow is linear in the vortex velocity and given by [17] 



2'Kha v, ; 



0. 



(2) 



Equation (2) can also be understood from the perspec- 
tive of momentum balance. A moving vortex imprints a 
phase discontinuity on the superfluid wave function. The 
Joscphson relation A/i = hdtA(p connects the resulting 
chemical potential to the time derivative of the relative 
phase difference, cf. Fig. 2(b). The chemical potential 
drop will be balanced by a flow of particles, which results 
in momentum transfer from the particles to the moving 
vortex, perpendicular to the vortex velocity v„ and pro- 
portional to its magnitude. The proportionality factor 
a relates the change in the system's momentum to the 
vortex velocity. 



The Hall conductivity can be related to the drift ve- 
locity of a vortex. From (1) and (2) we get 

v» = ^v s . (3) 

In a system with a low density of vortices we can ne- 
glect the effects of vortex-vortex interactions. Consid- 
ering strictly dissipation-less flow, we obtain from the 
Joscphson relation a semi-classical expression for the Hall 
conductivity 



q n s v s 
2iThn l ,v l , 



h n v 



(4) 



where n v is the density of vortices and q the boson charge. 

In systems with Galilean invariance, in a reference 
frame moving at the vortex velocity, there should be no 
forces acting on the vortex. This requires v„ = v s and 
therefore sets a = n s [18, 19]. This relation is modified 
in the presence of a lattice, as we discuss now. 

We consider the standard model for interacting bosons 
on a lattice [14] 



H 



t 



<r,r'> 



b\b,e iA * 



blh,e 



-iA r 



+ -1)- (5) 

r r 

where b\ creates a boson on site r, t is the hopping am- 
plitude, U the on-site repulsion; A rr r = q A • dx is 
the phase factor due to applied gauge field A. We work 
in units where h = c = 1, and likewise we set the lattice 
constant a = 1. 

We first note that vortices live on the center of the pla- 
quettes of the lattice. We explicitly construct operators 
T x and T y which translate a vortex by one lattice con- 
stant in the x and y directions. We show that they obey 
the commutation relation (see supplementary materials 
for a full derivation) 



2iriN b /N 



(6) 



Here, Nb is the total boson number operator, and TV is 
the number of sites. We denote the particle filling by 
rib = Nb/N. From Eq. (6), we see that the Berry phase 
acquired by moving a vortex around a dual lattice pla- 
quette is 



2ira = 2ir(nb + p) 



with 



pez. 



(7) 



The integer p arises from the 2-7T ambiguity in Eq. (6). 

Equation (7) can also be understood in terms of mo- 
mentum balance [20]. Following Paramekanti and Vish- 
wanath [21], we note that Eq. (6) implies that when a vor- 
tex is transported by Ay sites along e y , the momentum 
of the system changes by AP X = 27m;,Ay. At the same 
time we can integrate Eq. (2) to obtain AP X = 2TraAy. 
Combining the two results and taking into account that 
momentum is only conserved up to a reciprocal lattice 
vector 27rp leads us to Eq. (7). 
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The above results, together Eq. (4), imply a similar 
relation for the Hall conductivity, 

2 

(Txy n v = \ (n b + p) with p e Z. (8) 

While Eqs. (4) and (8) are a semi-classical derivation 
of the Hall conductivity, in Sec. IV we derive an exact 
relation between a and the Hall conductivity for a system 
containing one vortex, 

o- xy = ^Na = ^N(n b + p) with p e Z. (9) 

In the remainder of this paper we investigate the re- 
lations (7-9) throughout the phase diagram of the Bose- 
Hubbard model. In Sec. II we study these relations in 
the Gross-Pitaevskii and Mott transition limits. In Sec- 
tions III, IV we study the transition between different 
p-sectors in the hard core boson limit. We complete the 
phase diagram using numerical calculations in Sec. V. 

II. LOW-ENERGY LIMITS 

We start by discussing low energy limits of Eq. (5) 
where a diverging length-scale enables the derivation of 
a continuum low-energy theory. In these limits a and 
a xy can be deduced directly We review the derivation 
of the low-energy theories for weak (U t) and strong 
(U ^s> t) interactions. In both cases we start by rewriting 
(5) as a coherent state path integral with the following 
action for the complex valued field ipi 

s = f dT \ E^(^-^- t E(^^ e "^ 3+cx -) 
+ ¥E^i 2 (i^i 2 - 1 )}- ( 10 ) 

The Gross-Pitaevskii limit — In the weakly inter- 
act ing lim it, the Gross-Pitaevskii healing length £ GPE = 
ayt/Urvb is much larger than the lattice spacing a. This 
enables a straight-forward gradient expansion of (10). To 
lowest order in gradients we obtain the continuum action 

S = JdTdx{iP*d T 4> + a 2 t\(V -iqA)ip\ 2 + •••}• (11) 

Using the above expression we can now derive the co- 
efficient a in the Gross-Pitaevskii limit. When written in 
terms of i\) = ^Jn~b exp(ii^), the term ip*d T il) in Eq. (11) 
leads to a purely imaginary C T = liximrib, where m is a 
field that counts the winding of the phase $ [17]. Con- 
sider the action of a field configuration associated with 
taking a vortex around a closed loop of area S. A lit- 
tle reflection shows that the regions outside the loop do 
not change the value of m while those inside contribute 
unity per particle. Hence, C T gives rise to a Berry phase 
of 2nintS [18]. This observations fixes 

a = rib. (12) 



Around the Mott insulator At strong interac- 
tions and integer filling, rib G N, the Hamiltonian (5) 
stabilizes a localized Mott insulating phase with vanish- 
ing superfluid fraction ^ = (ipi) = [14]. In the insu- 
lating phase all sites are occupied by exactly rib bosons. 
Both the addition or removal of a particle is protected 
by a finite gap. This gap closes at the boundary of the 
Mott lobes in the phase diagram of Fig. 1, whereby at 
the lower (upper) boundary of the Mott lobe the hole 
(particle) gap vanishes. Hence, the tip of the Mott lobe 
represents a multi-critical point where the particle and 
hole gap close simultaneously [22]. In the following we 
focus on this multi-critical point. 

When both the particle and hole gap vanish, an en- 
hanced symmetry in the low energy sector emerges. In- 
stead of going through the standard procedure of deriving 
the low-energy theory from microscopic considerations 
[23, 24] we motivate the effective action via its symmetry 
properties. We expect the following particle-hole sym- 
metry (PHS) to hold tp and A -» -A [25]. To 
leading order in powers of ip we find 

S = JdTd X {sk- b \d T 4>\ 2 + a 2 tn b \(W - iqA)$\ 2 + ...}. 

(13) 

The gradient expansion leading to an effective continuum 
theory is controlled by the diverging correlation length 
close to the second order phase transition into the Mott 
insulating state. 

A direct consequence of PHS in the continuum theory 
(13) is <j xy (A) = <J X y{— A). Together with the Onsagcr 
relation a xy (A) = —a xy (—A) we obtain o~ xy = 0. This 
result can also be understood in terms of vortex motion. 
As opposed to the Gross-Pitaevskii action, the continuum 
theory (13) is real. Hence it does not give rise to any 
Berry phase when a vortex is moved around a closed loop 
and we conclude 

a = 0. (14) 

Starting from the PHS points, we expect to find lines 
with a = in the n/U—t/U phase diagram of the Bose 
Hubbard model. From Eq. (7) on the other hand, we 
know that at a fixed density, a can only change by an 
integer. This leads to the conclusion that the lines with 
a = arc bound to lines of integer fillings in the phase 
diagram, c.f. Fig. 1. 



III. HARD CORE BOSONS LIMIT 

We now consider the limit t/U and fx/U —> m, 
where m is an integer. In Fig. 1, these limits lie in- 
between two Mott lobes. The two states with m and 
in + 1 bosons per site are degenerate single-site states of 
Hamiltonian (5). States with different fillings are sepa- 
rated by a gap of order U, and do not appear in the low 
energy theory. 
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We use a Schrieffer-Wolf transformation [26] to project 
the Hamiltonian (5) onto the subspace with only m and 
m + 1 bosons per site. The resulting Hamiltonian cor- 
responds to hard-core bosons (HCB) and can be written 
using spin-i operators; Sf + \ is the onsite number op- 
erator, and Sf (S~) raises (lowers) the occupation from 
mtom+1 (and vice versa). At zeroth order in t/U , the 
HCB Hamiltonian is given by 

n$ = -(m + 1)* ]T (e iA "'S+S;, + h.c.) - ^ s z r . 

<r,r') r 

(15) 

The HCB Hamiltonian (15) has an emergent charge 
conjugation symmetry. One defines the unitary transfor- 
mation 



C = exp (i-K ^2 S$ 



(16) 



C transforms particles into holes, i.e., S Z C = —S*, 
and 



qA, -fi) 



(17) 



At half filling for the hard core bosons, the Hamilto- 
nian (15) is independent of \i and hence Eq. (17) implies 
invariance under A — > —A. Hence, the Onsager relation 
a xy (A) = —a X y(—A.) implies that for half integer fillings 
(rib = \ + m) 



0. 



(18) 



Note that the situation at the HCB limits and at the 
tip of the Mott lobes are qualitatively different. The Hall 
conductivity at the tip of the Mott lobe vanishes due to a 
zero-crossing of a when p = —rib- In the HCB limit, the 
integer p jumps exactly at rib = m+^. In other words: in 
the first case it reflects a particle hole-symmetry between 
rib — 1 and rib+1 while in the latter the symmetry connects 
rib — 1 and rib at rib = m+ ^. The symmetry at the HCB 
limits has a remarkable consequence for a xy in the full 
phase diagram of the model, as we shall now show. 



IV. AWAY FROM THE HARD CORE BOSON 
LIMIT 

We now consider the effect of a finite but small value 
of t/U. Second order processes in which a virtual excita- 
tion with an on-site occupation of m — 1 or m + 2 bosons 
arc created lead to corrections to the Hamiltonian (15) 
of order t 2 /U. Taking into account all the different pro- 
cesses we obtain up to irrelevant re-normalizations of the 
parameters in "H H c 



a/ 1 ) —i/v 

T-HC — rI -H 



(o) 



» E 

«r,r';r"» 



■'S+S^iSpr+^+h.C.), 

(19) 

where ((r, r';r")) denote sites r and r' which are nearest 
neighbors of site r", and e m = (m + l)(m + 2)t 2 /U. 



The new terms in Eq. (19) break the charge conjuga- 
tion symmetry, and therefore the Hall conductivity at 
exactly half integer filling does not vanish; below we cal- 
culate it in the limit of small t/U. 

We consider the model Eq. (19) on a torus of size 
N = L x L y , with N even. The gauge field A describes 
a uniform flux penetrating the surface of the torus. We 
take the total flux to be one flux quantum, which in- 
duces one vortex into the system. An important gauge 
invariant quantity described by the gauge field are the 
two Wilson line functions [27] 



■My) = j> dxA x , 



= j> dyA y 



(20) 



We define 6^ = $ x (y = 0) and Q y = $ y (x = 0). Chang- 
ing the values of Q x and Oj, corresponds to threading 
Aharonov-Bohm (AB) fluxes through the two holes of 
the torus [27]. 

The Hall conductivity at zero temperature, for a gen- 
eral many-body Hamiltonian can be calculated by inte- 
grating the Berry curvature [28] 



T is given by 



a* o (0) 



(21) 



(22) 



Here 4'o(®) is the many-body ground state wave function 
which depends on the Aharonov-Bohm fluxes through the 
holes of the torus. 

Remarkably, the Hall conductivity in the presence of 
one vortex can be calculated analytically at half filling. 
The key ingredients are degeneracies in the spectrum that 
occur for t/U = [27, 29] and serve as point (monopole) 
sources for the Berry curvature T [30]. 

To understand the degeneracies, we consider an effec- 
tive Hamiltonian for the vortex hopping between dual 
lattice sites. As shown in Ref. [27, 29], this is given by 



iA 



RRl btb^, + h 



R R' 



<R.R') 



£>(R-Rv)&k*R, (23) 



R 

where creates a vortex on a dual lattice site, and 
ty w t. The dual gauge field's flux is given by the boson 
density, VxA D = $ c = 2nn b . The potential [/(R-R v ) 
for the vortex position arises due to the fact that the Wil- 
son lines (20) break translational symmetry on the torus. 
In fact, as shown in [27, 29] for one flux quantum pene- 
trating the surface of the torus, all translational symme- 
tries are absent, and the potential U(r — Ry) acquires its 
minimum at a point Ry for which the Wilson lines both 
take on the value tt. 




FIG. 3. Vortex Hamiltonian. (a) The vortex hops betwe 
dual lattice sites (red) , where at half filling for the bosons 1 
average dual flux per plaquette is | a flux quantum. T 
circular contours represent the equipotential contours for i 
confining potential for the vortex. When the Aharonov-Bol 
fluxes are tuned to 0° , 0$) , the minimum of the potent 
is situated on a (direct) lattice site and the vortex grou 

state has two degenerate states \R), \L). (b) The two grou. 

states correspond to two different charge density wave con- 
figurations, with | R) (\L)) corresponding to depleted (excess) 
charge at Ry. When the particle hole symmetry breaking 
terms of Eq. (19) are introduced, hopping terms through Rv 
(depicted by the blue arrow) lower the energy of the state \L) 
relative to \R). 



If the point Ry lies on a site of the direct lattice, the 
cigenstatcs of Hy (in a symmetric gauge) can be writ- 
ten as ip(R- R v ) = /(|R- R v |)e 4m¥,( k~ Rv) . Here 
ip (R — Ry) denotes the angle between R — Ry and the 
x-axis, and m = 0, ±1,±2. At half filling, the average 
dual flux per plaquette is $£> = 7r and the ground state 
is doubly degenerate with m = 0,1. The two states \R) 
(m = 0) and \L) (m = 1) represent states with clock- 
wise and counter-clockwise vortex currents, respectively, 
as depicted in Fig. 3(b). Note that this two- fold degen- 
eracy occurs for N distinct values of 0. 

To calculate a xy we need to analyze the spectrum 
around the N degeneracy points. Around these points, 
the Hamiltonian restricted to the \R) and \L) basis is of 
the form Hy = h ■ a . To find h, we first notice that if 
©° = (0°, 0°) is a degeneracy point, tuning away from 
it by = 0° + A0, moves R v as [27] 

R v = Rv + AR v , Ai* v = -^^L a AQp, (24) 

where a, (3 = x,y and indices are not summed. Thus, 
tuning away from 0° breaks the degeneracy between the 
two ground states \R) and |L), as it shifts the minimum 
of the potential U(R — Ry)- Second, the terms cx e m 
in (19) lift the degeneracy even at = 0°; the assisted 
hopping terms through Ry (blue arrow in Fig. 3) favor 
\L) over \R). 

Together, these two effects give rise to the following 
low energy Hamiltonian for each degeneracy point (see 
supplementary materials for details), 

Hy = U (-AQ y a x + AQ x a y ) + la z . (25) 




FIG. 4. Berry monopoles. The parameter space of 
<d x ,G y , and e m . The yellow surface denotes the particle- hole 
symmetric manifold e m = 0, on which the red dots denote the 
degeneracy points for N values of (Q x , ©»)■ The blue surfaces 
denote e m < and e m > 0. The integral of the curvature T 
on these surfaces counts the number of sources in the three 
dimensional parameter space which they enclose. 



where the energy scales appearing above are U cx , 
and e cx e m . 

We use Eqs. (21-25) to calculate the Hall conductivity 
for one vortex. Consider the Hamiltonian Tile of Eq. (19) 
where we let the parameter e m take on both negative 
and positive values. For < e rn /t -C 1, the many body 
ground state ^ + (0,e m ) is non-degenerate, and likewise 
^-(0, e m ) for —1 <C e m /t < 0. For e TO = the two states 
become degenerate at a set of N values of space. The 
Berry connection of the ground state manifold at e m = 
is given by 




and must satisfy / d 2 T — due to particle hole sym- 
metry at e m = 0. As a result, 

>0) + < 0) = (27) 

Next, we consider the space of Q x ,Q y ,e m , which has 
the topology of a thick torus, as depicted in Fig. 4. We 
are interested in the integral of the Berry curvature T 
on the surfaces with e rn > and e m < 0, which yield 
o'xyi^m > 0) and a xy (e m < 0) respectively. At the same 
time, an analog of Gauss' law for T implies that the inte- 
gral of J- on these surfaces counts the number of sources 
for T [30, 31]. These are just the degeneracy points dis- 
cussed above, which are all described by Eq. (25) and 
therefore correspond to sources with charge +J. This 
leads to 

o- xy (e m > 0) - <J xy {e m < 0) = N. (28) 
Combining Eqs. (27) and (28) gives 

a xy {0 < t/U< 1) = f . (29) 
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Before concluding this section, we note that Eq. (24) 
leads to an exact relation between a and the Hall con- 
ductivity of one vortex. From Eq. (24), the Berry phase 
for moving a vortex around a plaquette is given by 



t/U = 0.2 



2-KOL 



d & T 



S(C) 



(30) 



where the contour C defines a square of size 2n/L x x 
2ir/L y in flux space, and S(C) is the surface it bounds. 
Therefore, for one vortex, 



aN = i- 



(31) 



Consider again the phase diagram of the Bose Hubbard 
model. From the above discussion, we conclude that the 
transitions lines between two integer values for p emanate 
from the HCB points, and move to higher densities with 
increasing t/U . These lines all correspond to changes of 
the integer p by unity. The PHS lines emanating from 
the neighboring Mott lobe tips terminate at the transi- 
tion lines, cf. Fig. 1. Together, they define regions with 
negative a and Hall conductivity. 



V. EVOLUTION OF THE TRANSITION LINES 

We numerically calculate the Chern number (21) for 
one vortex, to obtain the behavior of the integer p in 
the full parameter regime of the Bose Hubbard model. 
We use a Lanczos algorithm [32] to find the ground- 
state wave function \I>o(®) for different AB fluxes. Using 
a standard procedure [33] to numerically integrate the 
Berry curvature (22) we obtain the Hall conductivity for 
different values of t/U and n^. 

In Fig. 5 we show the results obtained for a 3 x 3 cluster, 
cf. supplementary materials. We indicate which integer 
p describes the Hall conductivity in panel (b). Panel (a) 
shows a trace of a xy for different particle numbers at 
t/U = 0.2. In both panels the region in the phase dia- 
gram where a xy < are marked with yellow hatches. As 
expected from the calculation at half-filling, the transi- 
tion lines between two integer value of p move to higher 
densities for t/U > 0. Remarkably, the transition lines 
intersect the integer density line at increasing values of 
t/U for higher densities. As a direct consequence the area 
of negative Hall conductivity increases for higher densi- 
ties. This is in contrast to the decreasing extent of the 
Mott insulating phases indicated by the yellow bars in 
Fig. 5(b). This surprising behavior of the Bosc-Hubbard 
model is explained below. 

In order to see at which values of t/U a sign change of 
a should be expected at integer fillings, we consider the 
healing le ngth ^ G PE , which sets the size of a vortex. For 
£gpe = a s/t/Unj, <C a the size of a vortex is much smaller 
then the lattice spacing a [34], and the Bose Hubbard 
model maps onto the quantum rotor model [10, 35] which 




FIG. 5. Numerical results, (a) The Chern number calcu- 
lated numerically on a 3 x 3 system at t/U = 0.2 for various 
particle numbers TV. The different branches are described by 
different integers p. The regions where the resulting Hall con- 
ductivity is negative is hatched in yellow, (b) Phase diagram 
obtained from cuts of the form of panel (a). Different colors 
indicate the different integers p. The lines where the integer 
p changes by one are the topological transitions. The Mott 
insulators at integer filling are indicated by the yellow bars. 



has an emergent PHS at integer filling [17, 36]. We have 
seen that PHS implies a = 0. The dependence of £ G pe 
on the mean-field interaction U therefore explains the 
growing extent of the negative Hall conductivity. This 
has to be contrasted to the bosonic enhancement factors 
in the hopping terms oc y/fib which facilitate the melting 
of the Mott insulator and lead to smaller Mott lobes at 
increasing densities. 

Finally, in Fig. 1 we present the numerical results as a 
function of t/U and \x/U . To translate from the results 
at fixed density n& to a fixed chemical potential p, we use 
a standard mean-field approach [37]. 



VI. EXPERIMENTAL VERIFICATION 

Our results have a direct experimental signature in 
terms of the vortex flow velocity in a moving system of 
lattice bosons. If o~ xy is positive (negative) the vortices 
move with (against) the supcrfluid flow, cf. Eq. (3). 

In a cold-atoms setup, the direction and speed of the 
vortex flow can be measured with in-situ imaging tech- 
niques [38]. However, in the strongly interaction regime 
the vortex core is smaller then a lattice spacing a. Hence, 
to make the vortex visible in the density profile, the sys- 
tem parameters have to be ramped into the weakly in- 
teraction regime before imaging. 
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The sign change of the Hall conductivity can also be 
measured studying collective modes in a trap. When the 
atom cloud is displaced from the minimum of the har- 
monic trap the atoms start to oscillate in the trap [39]. 
A non-vanishing a xy induces a transverse force on this 
dipole mode leading to a rotation of the axis of oscilla- 
tion. Depending on the sign of the Hall conductivity the 
rotation is clock or counter-clockwise. 



VII. DISCUSSION AND OUTLOOK 



In this paper we focused on vortex dynamics for the 
Bose-Hubbard model. We mapped the sectors corre- 
sponding to different integer p which characterizes the 
Hall conductivity and vortex motion throughout the 
phase diagram. We found that close to the Mott insulat- 
ing phases the sign of a xy is reversed and vortices flow 
against the applied current. 

Our results are obtained neglecting vortex- vortex in- 
teractions or disorder. We note, however, that there are 
an infinite number of particle hole symmetric points in 
the zero temperature phase diagram: at the tip of every 
Mott lobe, and in between two adjacent Mott lobes. We 
saw that the latter necessarily slice the full phase dia- 
gram into an infinite number of different p-sectors. This 
underlying structure cannot be removed by the inclusion 
of vortex-vortex interactions or disorder. However, the 
transition lines are expected to change their exact loca- 
tion and to be smoothed out by these effects, as well as 
by finite temperature. 

Incidentally, reversal of the Hall conductivity have 
been repeatedly measured in many strongly correlated 
electronic materials, including high temperature su- 
perconductors, e.g. in Refs. [3-5]. These experiments 
are beyond the direct applicability of our model. An 
extension to treat these materials is an interesting future 
direction. As discussed above, a clean verification of our 
predictions is possible in systems of cold atoms. 
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FIG. 6. The gauge choice Eq. (A4). Arrows on links represent 
the values for A r y : Black arrows represent a value of B, blue 
arrows represent BL X , red represent Q x /L x , and magenta 
Qy/Ly. Open links stand for periodic boundary conditions. 



Appendix A: Vortex translation operators 

We derive the commutation relations for the vortex 
translation operators 



%T y = T y % cxp(2niN b /N) 



(Al) 



Our derivation of equations (Al) holds for any interaction 
strength U in the Bose Hubbard model (5). We consider 
the Hamiltonian of Eq. (5) on a torus with N = L x L y 
sites. The gauge field A rr > describes one flux quantum 
piercing the surface of the torus uniformly, whereby the 
flux per plaquette is given by 



B = 



2tt 



(A2) 



An important gauge invariant quantity described by the 
gauge field A are the Wilson loop functions 



x (y) = f dxA Xl ® y (x) = f 



dyA y 



(A3) 



We choose a continuous paramctrization of the gauge 
field, which yields a continuous family of Wilson line 
function. Our gauge choice is given by 



Al 



r.r+x 



yBL x 8 x r 



(A4) 



where x = 0, L x — 1, and y = 0, L y — 1. Our gauge 
choice is shown in Fig. 6. The two Wilson loop functions 
are given by 

xBLy + Q y , (A5) 
(A6) 



%(y) = -yBL x + 6 
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Note that the parameters Q x ,Q y € [0,27r] define a con- 
tinuous family of Hamiltonians, which are inequivalent 
under gauge transformation. 

Let t x and t y be lattice translation operators, i.e. 



t x b r t x b r 



tyKty = b r + y 



Translating the Hamiltonian by t x and t y leads to 
tlU[A\t a =U[)S a \ 



A 



(<*) 



A 



t Q (r),t Q (r')- 



(A8) 



The gauge invariant content of the new gauge fields A( Q ) 
is the same flux per plaquette B as for A, however the 
Wilson line functions of Eq. (A3) arc shifted by one lat- 
tice constant, 



t x : <f> y (x) -> %{x) = <!> y {x + 1) 
t y : $ x {y) -> $ x (y) = $ x (y + i). 



(A9) 



Note that the values of the Wilson lines cannot be 
changed by a gauge transformation. Therefore if we con- 
jugate the Hamiltonian by t x or t y we cannot make a 
gauge transformation back to the original Hamiltonian. 
However, we can find gauge transformations U x and U y 
that yield the following relations: 



Uttin [A(0)] t a ui=n 

where = (<d x , Q y ), and 



A - A0 



(a) 



(A6^)^ = e a0 BL 



(A10) 



(All) 



where e a ^ is an antisymmetric tensor with e xy = 1. 
We parametrize the unitaries U x and U y as 

U x = cxp(i27r^Xr b t b r), U y = cxp(i2n X y r b\b r ) , 

r r 

where the functions x" are given by 

Xr = J dv' ■ (A(0-A0 (a) )-A (Q >) . (A12) 

We now calculate the commutation relation between 
T x and T y . An explicit formula for Xr using the gauge 
choice Eq. (A4) reads 

X x (r) = -BL x yS xfi , 

X v (r) = Bx. (A13) 

Multiplying these operators we get 

T y % = t y t x exp ^2 (Xr + Xr-x) b% j , 

%T y = t x t y CXp ^ (Xr-y + Xr) &r&r J 

= 7^7;exp(iT). (A14) 



In the above equation, we have used [t x ,t y ] = 0. The 
factor exp(iT) is given by 



T = uj r 6 r 6 r , 



(A15) 



( A7 ) with 



^r=Xr -Xr-y+Xr-x-X^- (A16) 



Using Eq. (A13) we get 



LU r = —B, 



(A17) 



and substituting this into Eq. (A14) we arrive at our final 
result 



%T y = T y % exp(2niN b /N). 



(A18) 



Here, Nf, is the boson number operator and N are the 
number of sites. 

Note that although explicit gauge choice were made in 
the derivation of Eq. (Al),the result is gauge invariant: a 
different gauge choice in Eq. (Al) yield the same result. 

To relate T x and T y to the vortex position, we note that 
the vortex position can only depend on the values the 
Wilson lines & x (y) and $ y (x), as these the only gauge 
invariant quantities which break the translational sym- 
metry on the torus. Therefore, the ground states of the 
continuous family of Hamiltonians "H[A(©)] corresponds 
to many body states \&(0) with vortex positions continu- 
ously parametrized by as well. The vortex position in 
the many body states \&(0) has quantum fluctuations; 
the amplitudes however are centered around a point Ry 
which depends only on [27, 29]. As a result, the action 
of Tx and T v shifts the vortex position by one lattice site 
in the x and y direction accordingly. 

To relate Eq. (A18) to 2ira, the Berry phase acquired 
by moving a vortex around a dual lattice plaquette, we 
note that 

exp(i27ra) = exp(i <j> dQ^A^) 

= <*(© o )|[/4lW2£/i|*(0 )), (A19) 

where A^ = i(\&(©)|9e M *(0)), and the line inte- 
gral is taken around a plaquette in flux space of size 
{BL y , BL X ). In Eq. (A19), U, = U(@\@ i+1 ) are adi- 
abatic evolution operators from Z to l+1 , and 

1 = 0° + A& {x) 

2 = 0° + A0 ( ' T) + A0 (y) 

3 = 0° + A0 fe) 



J 



0°. 



(A20) 

Using U 3 = TyU^Tj and U* = TjU%% and assuming 
that the ground state |*(0)) is non-degenerate, we find 

exp(^de^) = (*(0°)|rjr;v;i*(0 )) 

= cxp(-i27m;,) (A21) 
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FIG. 7. Vortex Hamiltonian. The green curves show the 
energy for a particle on a ring as a function of the flux through 
its center. At slightly higher (lower) flux than $_d = 7r the 
state m = 1 (m = 0) is the ground state. Two charge density 
waves centered at Rv lead to $u = ir+5, n — S and correspond 
to the states \L), \R) which have the same energy. Therefore, 
the state \L) (m = 1) corresponds to a charge density wave 
with excess density at the center site. 



Finally, from Eq. (A21) we see that the Berry flux 
through an elementary plaqucttc in flux space (which has 
the topology of a torus) is 27r(nb +p). All of the N el- 
ementary plaqucttc on the flux torus are identical. The 
Hall conductivity is given by the integral of the Berry 
curvature on the whole flux torus [28], and therefore, in 
the presence of one vortex, 

a xy = N(n b +p). (A22) 



Appendix B: Vortex hopping Hamiltonian 

In the following wc derive the form of Hamiltonian 
(25). Let us start with the terms arising from a change in 
the AB fluxes = 0° + A0, which moves the vortex po- 
tential J7(R — Ry) minimum according to (24). We can 
now apply degenerate perturbation theory in the sub- 
space of \R) and \L). The effect of A0 only leads to off 
diagonal matrix elements between \R) and \L) since both 
states have the same \ip(R. — Rv)| and the perturbation 
is diagonal with respect to R. The state \R) + e~ l{p \L) 
has excess weight along arctan(x/y) = ip, and therefore 
becomes the ground state for ARy along that direction, 



Hv = U (-AO v a x + Ae x cr v 



(Bl) 



We now consider the effect of the particle-hole sym- 
metry breaking terms of Eq. (19). We expect the two 
ground states \R) and \L) to conform to two charge den- 
sity wave orders centered at Rv, as the moving vortex 
exerts a force on the particles due to the Josephson re- 
lation. We note that the two charge density wave orders 



decay exponentially with the distance from Ry (the de- 
cay length scale is the lattice constant) [27, 29]. To see 
which of the states ( | R) or | L) ) has an excess (reduced) 
density at Rv, we consider an analogy with a particle 
hopping on a ring around Ry- 

In Fig. 7 we show the energy of a particle on a ring 
for the two states \R) with m = and \L) with m = 1 
as a function of the flux $o through the center of the 
ring. Note that the energy of \R) for = it — 8 is 
equal to that of \L) for <& D = n + 8 (at <J>£> = n the 
two states are degenerate). Consider now the vortex 
Hamiltonian of Eq. (23). At half filling for HCBs, we 
can consider two dual flux configurations which at Rv 
have $d(R v ) = n ± 8. Via $ d (t) = 2tt{S* + ±) they 
are related to two corresponding charge configurations 
(the two configurations are related by charge conjuga- 
tion). From the analogy to the particle on a ring, we 
can infer that these two charge configurations lead to the 
vortex ground states \L), \R) respectively. Wc therefore 
conclude that the state \L) (\R)) has an excess (reduced) 
density at Rv, respectively, cf. Fig. 7. 

While for the particle-hole symmetric point these 
charge configurations are equivalent energetically, the as- 
sisted hopping (t 2 /U) terms in Eq. (19) give different en- 
ergies for the two configurations. Due to the exponential 
decay of the charge density wave order [27, 29], the dif- 
ference between the expectation value of these terms in 
the states \L), \R) is also going to decay exponentially 
with the distance to Rv- 

To account for their effect wc estimate the energy 
change using mean field HCBs states for \L), \R) of the 
form |*) =n r K(tf r )), with 



|V0?r)> = cos(0 r )| |) +sin(0 P )e iv *| f>. 



(B2) 



Here ip r is the phase arising due to the votex. Due to 
the exponential decay of the charge density wave order 
we only consider the 3x3 cluster shown in Fig. 7. We 
evaluate the assisted hopping (19) in a state where the 
parameters i? r are chosen such that the center site has 
rib = 5 i 8, its nearest neighbors = \ =F 8/4, and 
the sites at the corners of the cluster have m = h. We 
find that the energy difference is dominated by assisted 
hopping terms which hop over the central site Rv 



E(\L))-E(\R))^-e m J2 



(L|e iA "'5+(^ v + i)5 r -;|L) 



-(R\e iA -'S+(S^+i)S-,\R) 



(B3) 



where r and r' are nearest neighbors of Rv • As discussed 
above, the state \L) corresponds to higher density at Ry, 
therefore the quantity above is negative. 

The combined effect of moving the vortex position Rv 
away from a direct lattice site and the PHS symmetry 
breaking terms of (19) leads to the low energy effective 
Hamiltonian near the degeneracy point 



Hy = U (-AQ y a x + AO x (r y ) + la z 



(B4) 
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which is given in Eq. (25). 



Appendix C: Exact diagonalization 

Wc calculate the ground state wave function for differ- 
ent AB fluxes using the ALPS Lanczos application [32] 
on a 3 x 3 cluster. We choose the same gauge choice a de- 
picted in Fig. 6. To obtain the phase diagram we truncate 



the local Hilbcrt-space to include all occupation states up 
to five particles per site. 

To get some insight into finite size effects we also cal- 
culate <7 xy for a 3 x 4 cluster at filling nt, = 1. To compare 
the two cluster sizes we estimate the Mott transition by 
considering the gap to the first excited state. We at- 
tribute the transition to a kink in the gap as a function 
ott/U. If we rescale the results for the Hall conductivity 
by the critical t/U the change from a xy = to a xy = 1 
obtained with the two clusters fall on top of each other. 
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